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ABSTRACT 

Surface brightness fluctuations have been proved to be a very powerful technique to determine the 
distance and characterize the stellar content in extragalactic systems. Nevertheless, before facing the 
problem of stellar content in distant galaxies, we need to calibrate the method onto nearby well-known 
systems. In this paper we analyze the properties at J and K s bands of a sample of 19 star clusters in 
the Large Magellanic Cloud (LMC), for which accurate near- infrared (NIR) resolved star photometry, 
and integrated photometry are available. For the same sample, we derive the SBF measurements 
in J and fsT s -bands. We use the multi-purpose stellar population code SPoT (Stellar Populations 
Tools) to simulate the color-magnitude diagram, stellar counts, integrated magnitudes, colors, and 
surface brightness fluctuations of each cluster. The present procedure allows us to estimate the age 
and metallicity of the clusters in a consistent way, and provides a new calibration of the empirical 
s-parameter. We take advantage of the high sensitivity of NIR surface brightness fluctuations to 
thermally pulsing asymptotic (TP-AGB) stars to test different mass- loss rates affecting the evolution 
of such stars. We argue that NIR-SBFs can contribute to the disentangling of the observable properties 
of TP-AGB stars, especially in galaxies, where a large number of these stars are present. 

Subject headings: stars: AGB and post-AGB — stars: carbon — stars: mass loss — galaxies: stellar 
content — (galaxies:) Magellanic Clouds — galaxies: star clusters 



1. INTRODUCTION 

In the study of galaxies, the age and chemical composi- 
tion of the stellar components are perhaps two of the ma- 
jor quantities to determine. Only then, we can trace the 
star formation, the chemical e nrichment and the assem- 
bly history of galaxies (see, e.g. lRenzi ni 2006). Since first 
applications, surface brightness fluctuations (SBFs) have 
been recognized to be effective to disentangle the evolu- 
tionary status of unre solved stellar populations in extra- 
galactic systems (e.g . iTonrv et alj|1990t fe lakcsl ee et al 



2001 



2003; 



iLiu et al]l2002l: ICantiello et al.ll200a Uensen et~ 



Cantiello et al 1 120051 ). SBFs are indeed much 



more sensitive to the brightest stars in the population 
in a given pas sband than integrated luminosities (e.g. 
ILiu et alJ feOOO), since they are defined as the ratio of the 
seco nd to the first moment of t he stellar luminosity func- 
tion (|Tonrv fc Schneidenll988l ). Thus, for instance near- 
infrared (NIR) SBFs may be efficiently used to detect the 
presence of intermediate-age stellar populations whose 
light is do minated by asymp totic giant branch (AGB) 
stars (e.g. iFrogel et al.lll990T ). Hence, it is now avail- 
able an additional tool to supply information on stellar 
systems other than the classical age/metallicity indica- 
tors based on integra ted light, as broadband colors, and 
spectral feature s ( e . g . I Wort he yj [T9 9 3 : Bres san et al.lll994l : 
lMarastonlfl998l ). 

Several works have been conducted to investigate the 
influence of AGB stars on integrated colors and spec- 
tral features in bot h metal-rich and metal - poor stel- 
lar p o pulations (e.g. Girardi fc Bertellil 119981: iMarastor] 
19981: (Brocato et all Il999at IMouhcine fc Lanconl 120021: 
Bruzual fc Charlotil2003t lMarastonll2005t iFagiolini et all 
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|2007| ). while only a few have been devoted to SBFs. At- 
tempts to understand how SBFs depend on the evolution 
parameters of AGB st ars have been m a de, in the opti- 
cal an d NIR bands, by ILiu et al.l (l2000f ) , iCantiello et al.l 
l|2003f) . and iRaimondo et all (l2005af) following differ- 
ent approaches. Liu and coworkers explored the ef- 
fects of changing, by an arbitrary factor, the lifetimes 
of post-Main Sequence (MS) stars in stellar popula tions 
older than 1 Gyr, whe reas ICantiello et alJ (|2003f ) and 
IRaimondo et all ((2005a) linked the lifetime of AGB and 
thermally pulsing AGB (TP-AGB) stars, in the age range 
0.1 ^ t(Gyr) < 14, to mass-loss process, the poor knowl- 
edge of which is among major uncertainties in modeling 
AGB stars. These stars may expel material at rates up 
to 10 -4 M yr _1 , eventually ejecting between 20 and 80 
per cent of their initial MS mass. Thus, their lifetime and 
the efficiency of the third dredge- up (TDU), which drives 
the formation of AGB Carbon-rich (C-type) star s, may 
be drastically reduced. IRaimondo et al.l (|2005al ) found 
that number variations of TP-AGB stars marginally af- 
fect optical SBF magnitudes for populations older than 
~ 1 Gyr, and, in general, less than ^0.5 mag in the 
NIR, confirming the reliability of the SBF metho d to 
measure distances to spheroids (jTonrv et al.l I200lh . In 
the case of young/intermediate-age populations, instead, 
SBFs appear to be highly dependent on the adopted 
mass-loss scenario, i.e. the number of TP-AGB stars 
(|Raimondo et a l. 2005a). An important caveat must be 
then understood when using SBFs to measure distances 
to galaxies suspected to host intermediate-age popula- 
tions. 

The high sensitivity of NIR-SBFs to the evolution- 
ary properties of AGB stars might be used, in turn, 
to test the prescriptions adopted in stellar models to 
describe the physical processes at work in such stars. 
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Big efforts to model stars beyond the core-helium burn- 
ing phase h ave been recently con ducted by several 
authors (e.g. iKitsikis fe WeisslfeOOTt iMarigo & Girardl 
l2007HCristallo et al.ll2008tlHerwigll2008r ). however, so far, 
a complete understanding of all aspects of the physics 
and numerical methods which properly describe AGB 
evolution is missing. As a matter of fact, the complex 
interplay between pulsation and mass-loss process, their 
dependence on metallicity and stellar mass, and inte- 
rior nucleosynthesis coupled with envelope chemical en- 
richment are the major sources of uncertainty in pre- 
dicting broad-band colors, abundances and luminosity of 
intermediate-age stellar populations. Therefore, what- 
ever the technique to derive the evolutionary status of 
the population is, consistency checks between model pre- 
dictions and observations for nearby resolved populations 
are still a prime step. 

In a previous paper (jRaimondo et al.ll2005ah we mea- 
sured optical-SBF amplitudes of 11 clusters in the Large 
Magellanic Cloud (LMC) using data from the Wide- 
Field Planetary Camera 2 (WFPC2) onboard Hubble 
Space Telescope (HST). In that work, we showed that 
the SBF data/model comparison suffers of a large uncer- 
tainty when a cluster contains a small number of stars 
evolved off the MS, especially giant stars. Thus, we sug- 
gested that a different algorithm has to be used to pre- 
dict SBF amplitudes for sta r clusters, due to the p aucity 
of stars in such systems. iGonzalez et al.1 ((2004) used 
a different approach: they built up 8 super-clusters by 
coadding clusters in the LMC and the Small Mag ellanic 
Cloud (SMC) having the same lSearle eTall (fl98fT SWB) 
class. The authors measured NIR-SBFs of these super- 
clusters using the Second Incremental and All Sky Data 
releases of the Two- Micron All Sky Survey (2MASS). 
To some extent, grouping clusters reduces stochastic ef- 
fects due to small numbers of stars on fast evolution- 
ary phases, however, this procedure suffers of the uncer- 
tainty in defining the membership of the cluster to the 
SWB classes as well as of the uncertainty of the rela- 
tionship b etween the s-paramete r and the cluster age. 
Moreov er, iMouhcine et al.l (|2005l ) and iRaimondo et al.l 
(|2005a[ ) ar gued that the discrepancy between predicted 
and observed SBFs of some super-clusters is likely due 
to observational reasons, as contamination by foreground 
sources, and low-quality photometry of stars in the cen- 
tral, most crowded, regions. Therefore, obtaining new 
SBF measures from NIR data having better spatial res- 
olution and sky stability is a necessary step toward con- 
straining more tightly stellar population models. 

In this paper we analyze the properties of a sam- 
ple of nearby populous intermediate-age stellar clusters 
expected to host TP-AGB stars. We s tudy 19 star 
cluste rs of LMC observed in the N IR by iFerraro et al.l 
(2004) and iMucciarelli et al.l (|2006f L and analyze their 
color-magnitude diagram (CMD), luminosity functions 
(LFs), integrated colors, magnitudes, and SBF ampli- 
tudes. From resolved-star photometry and total fluxes 
we measure SBFs in the J and K s pass-bands. The 
multi-purpose stellar popul ation synthesis code SP oT 
(Stellar Population Tools 2 , IRaimondo et al.l l2005aT ) is 
used to derive theoretical SBF amplitudes, CMDs, in- 
tegrated magnitudes, and colors of each cluster. We ex- 

2 http:/ /www. oa-teramo.inaf.it/SPoT 



plore the sensitivity of SBFs to TP-AGB stars by chang- 
ing the mass-loss rates. In the models presented here, the 
mass-loss rate is assumed to be the main parameter trig- 
gering the lifetime of TP-AGB stars, whereas the other 
relevant physical quantities quoted above are assumed to 
be similar in each models. The goal of the present paper 
is then twofold. First, we intend to reproduce the clus- 
ter's CMD, LFs, colors, and SBFs by means of a unique 
theoretical framework. Second, we address whether a 
connection between SBFs and TP-AGB star properties 
(number) and cluster age can be discerned in SBF data. 
The paper is organized as follows. § [2] details ingredi- 
ents of the stellar population synthesis (SPS) code, in- 
cluding a brief description of the method adopted to de- 
rive integrated quantities (colors and SBFs) . In § [3l we 
present the cluster sample and data, and in § [4] we com- 
pare synthetic CMDs, LFs, integrated colors, magnitudes 
and SBFs to the corresponding observed quantities. The 
final section provides a summary and discussion of our 
results. 

2. STELLAR POPULATION SYNTHESIS CODE 

The SPS models presented in this paper are based 
on an updated version of the stellar population synthe- 
sis code SPoT, designed to reproduce the properties of 
both resolved and unresolved stellar populations. SPS 
model predictions were applied and tested on star clus- 
ters of the Galaxy and MCs in previous papers (e.g . 
i Brocato et al.lll999al l2000l pOOl^aimondo et al.l l2002l. 
I2005al iCantiello et all 120031 120071 ). In this section, we 
briefly describe the main ingredients, techniques to de- 
rive integrated quantities, and outline changes. We refer 
the reader to the cited papers for more details. 

Differently from most SPS codes that are based on 
isochrones, the SPoT code directly relies on stellar evolu- 
tionary tracks. The presen t version utilizes the upd ated 
evolution stellar library by iPietrinferni et al.l (|2004T ) 3 for 
masses 0.5 < M/Mq < 11. The stellar models cover 
all the evolutionary phases from MS up to carbon igni- 
tion or the onset of thermal pulses (TPs). This allows 
us to analyze stellar populations in the age range from 
w 50 Mv r to 14 Gyr. S t ars of mass 0.1 < M/M© < 0.5 
are from IBrocato et al.l (|1998l ) . Masses lower than this 
limit (central H-ignition, M « 0.08 Mq) do not signifi- 
cantly contribute to the total mass of the cluster (see e.g. 
IChabrier fc Meral [T997). and so they are not taken into 
account. 

The mass of each star is randomly generated by using 
Monte Carlo techniques, while the mass distri bution is 
shape d by the Initial Mass Function (IMF) from lKroupal 
([20011 ) . The evolutionary line of each mass is then cal- 
culated by interpolating the available tracks in the mass 
grid. 

2.1. Horizontal Branch, AGB and TP-AGB stars 

The code is suited to simulate the color distribution 
of He-burning stars on the horizontal branch (HB) as 
a function of age, metallicity and mass-loss rate in old 
stellar populations (t <^ 5 Gyr; e.g. IBrocato et aflll999bl : 
IRaimondo et al.l 120021 ). The mass loss suffered by RGB 

3 The stellar evolutionary library by IPietrinferni et al. H2004T ) 
(BaSTI web-site: http://www.oa-teramo.inaf.it/BaSTil have been 
recently recomputed by the authors. 
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stars is evaluated according to the Reimers formulation, 
where the Reimers' parameter rj R (i refers to the i-th 
RGB star) follows a Gaussian distribution-function with 
a mean value (r) ) = 0.4 in our standard models. Thus, 
a generic RGB star losses mass at rate: 



M 



-4 ■ l(T 13 77f • LilU/Mi 



(1) 



where Li, Ri, Mi are respectively the star luminosity, ra- 
dius and total mass in solar units. The result of this pro- 
cedure is that He-burning stars are sprea d on the HB as 
observed in Galactic globula r clusters fe.g. lBrocato et all 
120001: iRaimondo et all [20021. SBF amplitudes of stellar 
populations older than ~ 5 Gyr were linked to t he HB 
morphology as a function of rj R in lCantiello et al. ( 2003). 

The evolution of stars off the zero-age HB till the 
double-shell pha se is followed by in t erpol ating the evolu- 
tionary tracks of lPietrinferni et alJ |2004) . In the case of 
intermediate-age populations, the stellar mass evolving 
off the turn-off (TO) point is high enough to prevent the 
RGB phase (i.e. no degenerate He core is developed), so 
that mass loss does not affect sizeably the color distribu- 
tion of He-burning stars. 

Beyond the early- AGB phase, stars with a MS mass in 
the range M ~ 1 — 8 M Q (the upper limit, called M up , 
depends on metallicity) undergo the TP- AGB phase . 
TP-AGB synthetic m odels (e.g. iRenzini fc VoI3 119811 : 
llben fc Renzinil 119831 ) and hybrid models, whi ch com- 
bine aspects of synthetic and full evolution (e.g. iMarigol 
1998), were developed with the aim to supply simpli- 
fied descriptions of stellar evolution in this phase, using 
recipes and descriptions based on full evolutionary mod- 
els and/or derived from empirical calibrations. In our 
code a specific routine evaluates the properties of ex- 
pecte d TP-AGB stars by integrating the a nalytic formu- 
las of W agenhuber fe Groenewegenl (11998J. here andaftcr 
WG98). as described in IRaimondo et al.1 (|2005al see Ap- 
pendix A). It calculates the time evolution of the stel- 
lar core mass and luminosity, the starting point being 
the helium-core mass, luminosity, and temperature at 
the first thermal pulse. The WG98 description includes 
three important effects: (i) the first pulses do not reach 
the full amplitude; (ii) the hot bottom burning process 
that occurs in massive stars; and (in) the TDU, even 
though this last phenomenon is treated in a simple way. 
In fact, in their models WG98 assumed that there is 
dredge-up only if the core mass is higher than a criti- 
cal value and that the free parameter A, describing the 
dredge-up efficiency, is constant. The limit of this for- 
mulation is that these descriptions and formulae do not 
account for the influence of abundance changes and nu- 
cleosynthesis in general on the evolution of the star. In 
the present work, the effe ctive temperat ure (T e ff) is eval- 
uated using relations of IWoodl (fl990l . hereafter W90), 
slightly shifted (AlogT e /y = —0.05) to match present 
evolutionary tracks. 

All stars are oxygen-rich (O-rich) when they enter the 
AGB phase. Whether or not they become O-rich stars 
primarily depends on the efficiency of the TDU occur- 
ring in the TP-AGB phase and o n the extent and time - 
variation of the mass loss (e.g. iStraniero et al.l [l997f ) . 
However, several effects play a role. In low-metallicity 
stars (Z < 0.004), the amount of oxygen in the envelope 
is so low that a few thermal pulses are sufficient to con- 



vert an O-rich star into a C-star (|Renzini fc Volilll98ll) . 
Moreover, the lower the metallicity the lower the min- 
imum mass for the onset of TDU (e.g. IStranier o et al.l 
2003). On the other hand, strong mass- loss episodes 
in TP-AGB stars may delay or even prevent the TDU 
occurrence and the formation of C-rich low-mass stars 
(jMarigo et al.l I1999T ). To introduce the metallicity de- 
pendence of the O- to C-star ratio, we consider the re- 
lat ion empirically derive d in galaxies of the Local Group 
by IGroenewegenl ([2007) . So that, we assume that the 
ratio of the duration of the two sub-phases for each star 
with a given metallicity is constant. 

Finally, post- AGB evolution experienced by stars be- 
fore entering the white dwarf cooling sequence has a neg- 
ligible impact in simulating the NIR properties of indi- 
vidual clusters, thus in the present work, such kind of 
stars are not considered. 

2.2. Mass loss of TP-AGB stars 

Major uncertainties in modeling the TP-AGB phase 
are related to the duration of this phase and to vari- 
ations of stellar luminosity and chemistry. As already 
stated, the duration is highly triggered by various physi- 
cal mechanisms, and primarily driven by the efficiency of 
mass-loss processes. To account for current uncertainties 
in the mass-loss rate determinations, we consider four dif- 
ferent mass-loss prescriptions, and assume that the other 
relevant physical quantities are described by the WG98 
and W90 formulations. Our aim is to investigate if and 
how they produce sizable effects on the foreseen observa- 
tional quantities, SBFs in particular. Our final goal is to 
test the sensitivity of SBFs to the number of TP-AGB 
stars, and to investigate whether SBF predictions/data 
of well studied MCs star clusters may give suggestions 
abo ut the number of TP-AG B stars. 

In IRaimondo et all (I2005al) we adopted the mass-loss 
rate formulation of lBloeckerl (fl995) based on the dynam- 
ical theoretical i nvestiga t ion on the atmospheres of Mira- 
like variables by[Bowen (1983) an d valid for long period 
variable (LPV) stars with periods P > 100 days: 



M B95 = 4.83 • lO' 9 ^ M- 21 L 2 7 M R 



(2) 

where the index i is suppressed for simplicity; M is 
the MS mass, r] B is a scaling factor, equal to unity in 
IRaimondo et a l. (2005a), and M R is the Reimers mass- 
loss rate with rj = 1. 

The four models considered in the present paper are 
defined by the following mass-loss laws: 



- Model A: The mass-loss rate of IBloecke'il (|1995l Eq. 
[5]) with rj B = 0.01 and rj B = 0.05 is applied to stars 
with M > 1.2Mq, whereas for lower masses the im- 
proved description f or cool winds of tip-AGB stars by 
IWachter et all (|2002f ) is adopted: 



M 



W02 



-4.52-6.81 log(T e// /2600)+2.471og(L/10 4 L o )-1.951og(^ 

(3) 



— Model B: The original Bloecker's law with r] B = 1 is 
applied to the full range of masses. 

— Model C: The em pirical mass-loss rate derived by 
Ivan Loon et al.1 (|2005f ) for O-rich dust-enshrouded AGB 
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log M 



V'05 



ing values as high as w 8 x 10 _5 MQyr _1 . A few stars 
with higher mass-loss rates are O-rich stars belonging to 
the LMC field population (small open circles). Note that 
here we plot models with a single metallicity, while data 
refer to galaxy populations, consisting of a mixture of 
stars with different age, metallicity and mass. There are 
5.65+l.Q51og(i/10Q00i Q )-6.31og(T e/ //3500iQ ) stu dies that conclusively show t hat the mass-loss rat e 



stars in the LMC is assumed. Through the modeling 
of spectral energy distributions, the authors estimated 
a mass-loss rate as a function of stellar luminosity and 
effective temperature: 



(4) 

— Model D: The law proposed by Ivan Loot] (l200l for 
very bright stars with dust-driven mass-loss is adopted 



M 



V06 



= 1.5 x 1(T 9 Z[Z Q }-°- 5 L[L G 



10.75 ^0.75 



(5) 



Here, Ay is t he visual extinc tion ranging from about 
0.01 up to 100 (|van Loor]|2007h . In our synthesis models 
we assumed Ay as a free parameter, and find a pos- 
teriori that values in the range 0.01 — 1 are suitable to 
reproduce present NIR data. This is because most dust- 
enshrouded AGB stars are visible at longer wavelengths 
compared to the ones used here (J and K s bands). More- 
over, from the analy s is of a sample of oxygen-rich AGB 
stars, iHeras fc Honvl ((20051 ) determined a mass-loss rate 
in the range 5 x 10~ 8 -^10~ 5 M Q yr _1 and a visual optical 
depth of 0.03 -f- 0.6 for stars with a small amount of dust 
in their envelopes. 

To illustrate differences between models, we computed 
the mass-loss rate evolution of three TP- AGB stars of 
M = 1.4, 2.8, and 4.6 M at fixed metallicity (Z = 
0.008). The resulting evolutionary paths as a function 
of the pulsation period are reported in Figure [TJ Theo- 
retical fundamental periods were calculated from period- 
mass-luminos i ty rel ations for long p eriod variables of 
IFox fe Woodl (| 19821 ). as reported by iMarigo fc Girardil 
(|2007l ). under the assumption that stars pulsate every- 
where on the TP- AGB phase, even if an instabili t y stri p 
should be considered (|Groenewegen fe de Jond Il994f) . 
Theoretical sequences are qualitatively compared with 
empirical mass-loss rates, derived from IR fluxes and CO 
radio line emission, and pulsation periods from observed 
light curves of a sampl e of obscured C- and O-rich stars 
(IWhitelock et al.ll2003f ) and C-stars (jGroenewegen et all 
120071) in LMC. In th e figure Galactic AGB s t ars ob served 
bv IWhitelock et all (l2o5l . iWinters et all (|200l . and 
iSchdier fc O lofsson (200J]) are also plotted. As remarked 
bv IWhitelock et all (|2006l ) . the two populations do not 
show remarkably differences in color and period relations. 
A similar conclusion can be inferred for the mass-loss 
rate-period relation, since at a given period mass-loss 
rates for LMC stars are in agreement with those obs erved 
in Galactic stars (see also lGroenewegen et alj [2007). De- 
spite a large spread at shorter periods, data suggest that 
the mass-loss rate increases strongly with period greater 
than logP ~ 2.6. A less s teep dependence can be argue 
at lower periods (see also iStraniero et all 120061 ). even if 
data (Galactic stars mostly) are more spread, as they 
suffer of the uncertainty on the distance, which appears 
in the expression used to derive mass-loss rate (see e.g. 
IWhitelock eFaIll2006l ). 

From Figure Q] it appears that in the case of model A 
(panel (a)) the evolutionary lines nicely reproduce the 
observed trend. As a matter of fact, the bulk of stars are 
fairly reproduced by stars having masses approximatively 
in the range ~ 1.4 -j- 4.5A/ Q , and mass-loss rates reach- 



explicitly depends on metallicity (|Ziilstra et all [19961) . 
however metallicity may have effect on grain growth, 
the number density and size of grains, thus altering the 
mass-loss rates from the stars. Finally, mass-loss rates 
may strongly depend on stellar pulsation. Stars of high 
metallicity, e.g. Z = 0.02, and mass are able to reach 
lower effective temperatures and higher mass-loss rates. 
In our models the central star typically has T e ff in the 
range ~ 2500—3800 K, in fair agreem ent with estimations 
found bv iGroenewegen et all (|2007f ) and Ivan Loon et aT] 
(2005) for MCs field populations, even thought the lat- 
ter sample includes cooler stars. If r\ B = 1 (model B) 
the mass-loss process is very efficient in stripping away 
the stellar envelope, and the evolution results s teep and 
rapid. The formulation of lvan Loon et all (|2005l ) applied 
to our TP- AGB modeling produces a rapid increase of 
the mass-loss rate since the beginn ing of the phase (Fig- 
ure [TJ model C). The relation from Ivan Lo^((20CJ7T) pro- 
duces a less steep behavior and longer periods at a given 
mass-loss rate (Figure [TJ model D). However, we recall 
this formulation critically depends on Ay. In the rest of 
the paper, unless explicitly stated otherwise, we present 
results of models A and C, even though computations 
and comparison for models B and D have also been per- 
formed. 

As conclusive remark, we notice that changing the 
mass-loss rate prescription may also affect the TP- AGB 
star structure, so that here we have performed a quali- 
tative comparison between results of our synthetic mod- 
eling of TP- AGB stars and observations. The detailed 
treatment of the evolution of stars along the complex TP- 
AGB phase requires accurate TP-AGB evolution stel- 
lar models, for instance the modeling of the two stages 
(O- and C-rich phases) as a function of stellar mass and 
metallicity, by taking into account the internal structure 
evolution, stellar pulsation, mass-lo ss mechanism, and 
dust and grain format i on (see e.g.. Ilzzard et al.1 120041 : 
IMarigo fc Girardill2007t iCristallo et alJl2009f ). Neverthe- 
less, the agreement between the data and the general pic- 
ture that we have presented is already very encouraging 
to investigate how these scenarios affect SBF predictions, 
as we do in the following sections. 

2.3. Integrated magnitudes 

To compute the cluster total fluxes we assume that the 
integrated luminosity is dominated by light emitted by 
their stellar component. On this basis, the total fluxes 
mainly depend on two quantities: 1 ) the flux fi emitted 
by the i-th star of mass M, age t, luminosity L, effective 
temperature T e ff, and chemical composition (Y, Z): 



fi[L(M, t, Y, Z),T eff (M, t, Y, Z),Y, Z] 



(6) 



and 2) the number of stars of mass M in a population 
counting N stars, $(M,N). fi is defined by the stel- 
lar evolution library and the temperature-color transfor- 
mation tables, while $(M, N) is strictly related to the 
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IMF. The Monte Carlo procedure we adopted is essential 
to simulate poorly populated stellar systems, like those 
we study in the following. It ensures that fast, under- 
sampled evolutionary phases, are correctly treated from 
a stochastic point of view. The above quantities are com- 
bined to calculate the total integrated flux F in a given 
photometric band: 



N 

F[N,t,Y,Z] = Yj i . 



(7) 



To account for stochastic effects, we computed a large 
number (N S i m ) of independent simulations for each set of 
population parameters (age, chemical composition, clus- 
ter mass, etc.). So that, we obtained statistical distri- 
butions of magnitudes (colors), produced by stochastic 
variations in the number and properties of bright and 
rare stars. For the purposes of this paper, N sim = 200 is 
appropriat e to reproduce the obser ved clusters (see also 
Figure 2 in IRaimondo et al.H2005al) . 

The conversion from theoretical quantities to mag- 
nitudes and colors is based on the BaSeL library 
(|Westera e t al. 2002, and reference therein), with the ex- 
ception of C- a nd O-rich stars, whose colo rs are derived 
from spectra of lLancon fc Mouhcind (120021). Th e hom og- 
enized photometric system of iBessell fc Brettl (|1988| ) is 
adopted, so that in the following magnitudes and col- 
ors involving the K filt e r are transformed into K 8 filter 
according to lCarpenterl (|2001l ). This assumption is a de- 
quate as extensively explained in iPessev et alj ((2006). 

2.4. Surface Brightness fluctuations 

If the type and flux of stars that comprise the stel- 
lar population luminosit y function are known, th e SBF 
magnitude is defined as (|Tonrv fc Schneider 1988): 



M = -2.5 log [^f] 

Z^i=l Ji 



(8) 



where M relies upon the Poissonian statistics. This is 
the case of our synthetic stellar population models, in 
which the flux of each star is known. In the case of 
a large number of independent simulations, as adopted 
here, the average SBF magnitude can be defined as 



(9) 



and the statistical uncertainty is derived as the standard 
deviation of the Mj distribution (j = l,N S i m ). 

3. CLUSTER SAMPLE AND DATA 

Th e NIR photometry of 19 clusters in the LMC i s 
from iFerraro et all (|2004h and iMucciarelli et alj (|2006l ). 
obtained with the SOFI imager/spectrometer mounted 
on the ESO 3.5m NTT. Each cluster typically con- 
tains 1000-1200 objects, with the exception of a sub- 
sample of poorly populated clusters (such as NGC 2190, 
NGC 2209, NGC 2249, NGC 1651, NGC 2162, and 
NGC 2173) and an handful of rich clusters (such as 
NGC 1783, and NGC 1978) with more than 1900 stars. 

The clusters are listed in Table [1] to gether with some 
properties. The age estimations from IMucciarelli et all 



(2006) ha ve been derived thro ugh the s-parameter cali- 
bration of IGirardi et al.l (|1995f ). based on stellar models 
computed with a certain amount of convective overshoot. 
[Fe/ H] determinations from either spectroscopic or pho- 
tometric data are also reported. The last three columns 
list the number of C-stars i dentified by v a rious authors 
and the s-parameter from Elson fc Falll (|1985l ). The 
membership of C-rich giant stars is not well established in 
some cases. For instance, a giant star in NGC 2136/37 is 
identify as a C-star wit h K a 10.7 mag and J — K = 1.53 
mag since the work by lAaronson fc Mouldl (Il985f). even 
its membership is not confirmed (jFrogel et al.lll990l ). 

We fixed the total F-magnitu de of each cluster to b e 
equal to the value published bv lGoudfrooii et ail (|2006f) . 
who measured integrated-light photometry in Johnson- 
Cousins V, R and / for a sample of 28 star clusters in the 
MCs. From their Table Al we selected measurements at 
an aperture radius corresponding to 1'.5, the same area 
used to the CMD analysis (see below) . For c luster s not 
included in the sample of iGoudfrooii et al. (2006]), we 
consider measurements by Ivan den Berehl (119811) . bear- 
ing in mind that they generally refer to a smaller aperture 
size (~ 50-60"). After a rescaling, the formal difference 
between the two sets of measures is n egligible, being of 
the o rder of (Vg — Vv) = —0.03 mag (jGoudfrooii et al.l 
12006ft . 

The distance modulus of LMC is settled to be (m — 
M)q = 18.40 ± 0.10 mag, on the basis of recent esti- 
mations (e.g iGrocholski et alJ 120071 : iTesta et al.l 120071 : 
iWalker et al.ll2001f l. It is worth mentioning that a varia- 
tion of the LMC distance within the above uncertainty, 
e.g. the adoption of the usual distance value (m — M) = 
18.50 mag, does not affect the results presented below. 
Similar consideration can be do ne for the Galactic e xtinc- 
tion law, here we use Ry = 3.4 (|Gordon et al.l [2003). The 
conversion of Ay to extinction coefficient s in other pho- 
tomet ric bands was done using formulae in lCardelli et alj 
(fl989h . 

4. DATA/MODELS COMPARISON 

The theoretical framework described in § [2] was 
adopted to derive fully consistent synthetic CMDs and 
integrated quantities. The comparison with resolved fea- 
tures of well studied star clusters (§ [3]) provides the op- 
portunity to calibrate NIR models in detail, by fitting 
the observed CMDs, star counts and luminosities, and 
at the same time integrated and SBF magnitudes. We 
take advantage of this joint analysis to infer indication on 
the efficiency of SBFs to estimate ages and metallicities 
of unresolved intermediate-age systems and to constrain 
the number and photometric properties of bright red gi- 
ant stars. 

4.1. CMD, luminosity function and stellar counts 

The procedure to fit models to data can be summa- 
rized as follows. First, a series of synthetic CMDs were 
computed to estimate the age and chemical composition 
by fitting the observed CMD, and the cumulative lumi- 
nosity function (CLF) of each cluster. The age step is 
10% of the expected age, while the metallicity is given 
by the available grid of stellar evolution models, i.e. no- 
interpolation was done. Each synthetic CMD includes a 
fraction of detached binaries whose mass ratio was dis- 
tributed homogeneously between 0.7 and 1, as found in 
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several LMC clusters (e.g. NGC 1818: lElson et al l 



1998 



NGC 1866: iBrocato et all [20031 : iBarmina et all 12002 : 
NGC 2173: iBertelli et al.ll2003l ). The fraction of binaries 
is a free parameter ranging from 10% up to 70% of the 
total. Age and metallicity evaluations are given by simul- 
taneously fitting the distribution (CLF) and photometric 
properties (CMD) of stars in the cluster. Note that the 
computed cluster's ^-magnitude (My*) is equal to the 
observed value. For the hand ful of clusters whose M\f l 
comes from the compilation of lvan den Be rgh (198]]), the 
systematic difference of the aperture-size was taken into 
account. Once the best-fitting parameters (age, metal- 
licity, and fraction of binaries) were established, we com- 
puted a series of 200 synthetic CMDs in order to take 
into account stochastic effects (see § 12.31 and § I2.4|) . It is 
relevant to emphasize that random extractions of stellar 
masses are fully independent, even though the same set 
of input parameters is assumed. Moreover, when sim- 
ulating the CMD we did not draw the same number of 
stars as observed, instead the number of stars placed in 
different regions of the synthetic CMD turns out from 
computations when the total magnitude is assumed to 
be equal to the observed one. 

Figure [5] shows observed (left) and synthetic (right) 
[K s , J - K S ]-CMD of each cluster. The synthetic CMD 
plotted is one of the 200 simulations computed in the 
framework of model A. C-type AGB stars are identified as 
stars, while M-type stars as asterisks. The vertical dot- 
ted lines indicate the rough separation between C-rich 
(redder) and O-rich (bluer) stars on the base of their 
J - K s color: (J - K s ) w 1.4 mag at Z > 0.008 and 
(J — K s ) 1.3 mag a t Z = 0.004 , as empirically deter- 
mined bv lCioni et all (pOOll l2003h in the LMC and SMC 
and confirmed by the cross-corr elation with spectroscopi c 
surveys for example in SMC bv lRaimondo et a l. ( 2005b). 
The clusters are arranged for increasing age, from 50 Myr 
to 5 Gyr. Metallicity ranges from Z = 0.02 down to 
Z = 0.004. The reddening value (Eb-v) labeled in each 
panel is the mean value of the measures available in the 
literature. 

As pointed out by iMucciarelli et al.l ((2006) and 
iFerraro et all (|2004f ) a significant contribution of LMC 
field stars is present in all clusters (see their figures 1-4, 
and figures 2-4, respectively). However, in young clusters 
a blue sequence is clearly visible at — 0.3 ^ ( J—K s ) Hz 0.3 
and K s ^ 15.5, corresponding to the brightest end of the 
cluster MS. At K s ~ 14 mag the region occupied by He- 
burning stars is also recognized. So that, it is possible 
to confidently identify cluster's stars as those placed at 
the left of the solid lines reported in Figure [2] (first six 
plots), notwithstanding a residual uncertainty remains 
to establish the membership of stars with K < 12 — 13 
mag (likely AGB stars). In the older clusters the obser- 
vational limit (K s « 18.5) does not permit to detect the 
TO point and for ages greater than ks 400 Myr it is not 
possible to easily identify field stars, since they overlap 
the cluster population. For all these reasons, to mini- 
mize the field contamination and to better estimate the 
cluster age and metallicity, the fit-procedure was applied 
to stars as far as r — 1.5' off from the cluster center 
(la rge circles in the left panels of Figure [2]), as suggested 
by IFerraro et al.l ((2004) . In addition, we used the pho- 
tometry of the control field a few arcmin away from each 



cluster to subtract the contribution of field stars to the 
observed CLF, after to have normalized and rescaled the 
data to the selected cluster's area. The resulting CLFs 
are used in the fit-procedure. 

Figure [3] compares observed and synthetic CLFs. Syn- 
thetic CLFs well reproduce the observed counterparts, 
especially at K s ^17—17.5 mag. This is a reli- 
able limit free by crow ding effects and completeness 
(|Mucciarelli et al.l l2006l ). For two clusters (NGC 1806 
and NGC 1783) the agreement is limited to the bright- 
est part of the CLF (K s J; 14 — 15 mag). Recent opti- 
cal observations have suggested the presence of two stel- 
lar populations in the area of the se two clusters, wit h 
an age difference of « 300 Myr (jMackev et all 120081 ) . 
The NIR data used here do not permit to reach the TO 
point, however our best age estimation for NGC 1783 
is about 1.2 Gyr, in fair agreement w i th the value re- 
cently suggested by IMucciarelli et al.l ([2007, t ~ 1.4 
G yr) and slightl y youn ger than the interval suggested 
bv lMackev et all (|20(E 1.6-2.2 Gyr), both derived from 
isochrones which include some overshooting. Other few 
clusters are recognized to be complex systems, for in- 
stance NGC 2136 and NGC 1987. The for mer is a com- 
pone nt of a potential triple cluster system ijHilker et al.l 
1995). The latter is a very poorly-popu lated cluster lying 
in a highly field-co ntaminated region (jCorsi et al.l I1994T 
IFerraro et all l200l . an occurrence that appears to be 
confirmed in Figure [3J Corsi and co-workers in their op- 
tical CMD discriminated three different populations: a 
very old field population, a young population producing 
a blue plume extendi ng the bright MS, a nd the cluster 
population. Moreover. iFrogel et al. (1990) identified two 
C-stars in the cluster area, whose membership, however, 
is determined uncertain for one of the two. 

As our models A-B-C-D differ only for the treatment 
of the TP- AGB phase, there is a good agreement between 
observed and synthetic CMDs till the early-AGB either 
for model A or models B, C and D. Thus similar ages 
and metallicities are predicted from all the models. In 
contrast, the four models produced a different number of 
stars in the subsequent TP- AGB phase (see below). 

Table [2] summarizes the age and metallicity, the mean 
My and Mk s absolute magnitudes, and J — K s and 
H — K s colors {models A) for each cluster. As we used 
the available metallicity grid without making any inter- 
polation, one can consider as metallicity uncertainty the 
difference between the best-fit value and the nearest two 
values in the grid. 

From the age estimations in Table [2] and s values from 
Table [T] we obtain the following age-s parameter calibra- 
tion (Figure HJz): 



Xogt = (5.37 ± 0.12) + (0.099 ± 0.003) x s (10) 

Our calibration based on can onical stellar evoluti onary 
tracks is compared to that of iGirardi et al.l (|1995l ) and 
ages from the literature in Figure \4}>. The use of both 
CMD and CLF has permitted to derive age estimations 
in agreement with those from the literature, however a 
future extension to the TO region will help in better de- 
fine the age of some clusters, especially those for which 
recent deep optical-data have revealed the presence of 
multiple populations (e.g. NGC 1806 and NGC 1783, 
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iMackev et al.l [2008) . Differences between the two cali- 
brations are mostly due to the e fficiency of overshoot- 
ing: in stellar conne ctive H cores (|Barmina et al.l l2002t 
iBrocato et alJl2003l and references therein). 

We compare now the expected number of C-rich AGB 
stars derived from models A-B-C-D to the observed one. 
In Figure [5] the number of C-rich AGB stars normalized 
to the total visual light is plotted as a function of the TO 
mass (Mto, Figure[5ji) and the s-parameter (Figure^). 
The error budget of data accounts for the uncertainty of 
integrated visual light and the membership of C-stars, 
while theoretical uncertainty is estimated taking into ac- 
count statistical fluctuations in the number of TP- AGB 
stars and the uncertainty of integrated visual light. As 
shown in Table [U the number and luminosity of bright 
AGB stars may suffer a large statistical fluctuation in 
clusters of similar age. This is due to the small number 
statistics, as each cluster contains at most a few TP- AGB 
stars. This is visible in Figure(S]i at Mto ~ 1-8 — 2.2 M , 
where a few clusters having nearly the same age posses 
a different number of TP- AGB stars. Around these val- 
ues there is a maximum in the frequenc y of the TP- AGB 
stars (see also iMarigo fe Girardil 120071 ). The qualitative 
analysis of Figure [5] also shows that the predicted trend 
is similar in all models, even thought the observed max- 
imum is nicely reproduced by model A, and to a lesser 
extent model D and C . We find that in the mass range 
1.7 < M/M Q < 2.5 and for metallicity Z = 0.008 the 
TP-AGB lifetime is of the order of ~ 2-3 Myr and re- 
duces both at lower and higher masses. This result is 
in agreement with the analysis recently performed by 
IMarigo fc GhaTdl ((20071) . 

4.2. Integrated NIR colors 

TP-AGB stars affect NIR-bands integrated luminos- 
ity and its uncertainty in low- metallicity intermediate- 
age massive clu s ters (see e.g. IGirardi fc Bertellil 119981: 
Marastonl Il998t iMouhcine fc Lanconl 120021 : iMarastonl 
20051) . In particular, NIR colo rs may be dominated 
by a handful of red giant stars dSantos & Fro£e][[T997| 
IBrocato et aHll999atlc7rvino fc Valls-Gab aud 200|). At 
the typical age of Galactic globular clusters (GGC), the 
presence of C-rich stars becomes more uncertain, as in 
GGC AGB stars are all observed to be O-rich, so that 
carbon does not appear to have be en dredged up into 
the e nvelop during thermal pulses (jHabing fc Olofssonl 
2003). It is worth noting that the assumption on mass 
loss is expected to have a marginal effect on integrated 
magnitudes and colors of very faint populations, where 
bright TP- AGB stars are statistically less frequent, o r 
even absent ([Santos fc FrogeHll997tlFagiorini et a"Lll2007l) . 

The differences between observed and predicted mag- 
nitudes (colors) of each cluster are plotted for models A 
and C in Figure [6] and Figure respectively. Theoret- 
ical values are calculated as the average over 200 sim- 
ulations (see § I2.3p . and error-bars are related to the 
standard deviation (ler). Simu lated values of My 1 are 
co mpared to mea s ureme nts of iGoudfrooii et al.l (2006) 
or Ivan den Ber gh (1981) (left higher panel), while pre- 
dicted values of Mk s and colors (J — K s an d H — 
K s ) are compared to data from iPessey et al.l (120061. 
top). IMucciarelli et alj (|2006l middle). and lPersson et al.l 
(|1983l bottom). Observations are decontaminated for 
the field contribution, and dereddened by using the mean 



Eb-v values labelled in Figure O 

Figure [6] shows that the agreement is nice f or all 
the clusters in the case of IPessev et all ([2001 ) and 
IMucciarelli et all (120061) . while a discrepancy can be 
seen with IPersson et al.l (jl98 3) . As pointed out by 
IPessev et al.l (j2006l ). a major discr epancy between their 
color measurements with those of IPersson et al.l (|1983f ) 
are caused by centering problems and to the presence of 
several relatively bright stars in the background field. In 
the figure it is clearly shown that colors of low-luminosity 
clust e rs have a large dispersion (see alsolSantos fc Frogell 
[19971: IBrocato et all [l999at iFagiolini et al.ll2007h . This 
is the case of the faint cluster NGC 2209, for which 
Pessev et all (120061) recognize the major difference with 
Persson et al.l (1983). For this cluster we estimated an 
age logi = 8.95 ± 0.10, in agre ement with our previou s 
value from optical-HST data (jRaimondo et alJl2005aD. 
The c luster contains t wo C-stars, identifie d by [Walker 
(|1971l ) and confirmed bv lFrogel et al.l (|1990l ). dominating 
NIR fluxes. No M-giants are possessed by the cluster. 

Inspection in Figure [7] shows that predicted K s magni- 
tudes still agree with observations within the errorbars , 
esp ecially in the case of da ta from IPessev et al.l (2006) 
and IMucciarelli et~a l. (2006). With a few exceptions, dif- 
ferences between observed and predicted colors are con- 
tained in the region marked by the dotted lines, illus- 
trating all the color values covered by simulations. Thus, 
synthetic colors are still able to reproduce the data in the 
case of models C, even if the agreement worsens at redder 
colors (i.e. J — K s It 1 mag and H — K s lt 0.3 mag), as 
a consequence of the smaller number of TP-AGB stars 
predicted by model C compared to model A. Nearly the 
same results can be inferred from the analysis of colors 
from model D, and slightly worsen for model B (see Fig- 
ure [8]) , where for sake of brevity only differences between 
observed and predicted if s -magnitudes and (J — K s ) col- 
ors are reported. 

In conclusion, observed NIR magnitudes and colors of 
LMC clusters are reproduced by model A, and to a lesser 
extent by model C, D, and B (the worst case) , within er- 
rorbars or at least within the color distributions obtained 
from all the simulations (200) computed for each cluster. 
Therefore, even though NIR integrated magnitudes and 
colors are sensitive to TP-AGB stars, they do not appear 
to be efficient enough to discriminate among physical as- 
sumptions which affect the number of bright TP-AGB 
stars in clusters, at least for the degree of accuracy of 
models proposed here (§ 12. 2|) . 

4.3. Surface Brightness Fluctuations 

In this section we present observed SBF magnitudes 
of individual clusters and compare them with theoretical 
predictions. In the case of star clusters the SBF measure- 
ments are obtained from the second moment of the stellar 
luminosity f unction directly derived from resolved-star 
photometry (lAihar fc Tonrvl 11994 iGonzalez et al.ll2004l : 
iRaimondo et al.l l2005a| ). Let us recall that this proce- 
dure measures the same physical quantity obtained for 
galaxies, and derived from the pixel-to-pixcl brightness 
variation in the galaxy image. Studying SBFs of LMC 
star clusters is then a relevant step to understand the 
SBF signal from unresolved galaxies where light is (to- 
tally or partially) emitted by the underlying mixture of 
old and intermediate-age populations. 
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SBF magnitudes of clusters were estimated as fol- 
lows. We coupled individual stellar pho tometry and total 
fluxes of the clusters both provided by iMucciarelli et al] 
(2006). The denominator is the cluster total flux, free 
from field contamination and reddening. The depth of 
the stellar photometry used here, contributing to the nu- 
merator, i s adequate for our pu rposes. In fact, as dis- 
cussed by lAihar fc Tonryj (p.994), in the optical and at 
least for halo GGGs, the second moment of the lumi- 
nosity (numerator), converges quickly, with 99% of the 
sum being obtained with the three brightest magnitudes 
of cluster stars. This is even more evident in NIR bands 
((Gonzalez et al.|[2004f) , where the total light is dominated 
typically by RGB and AGB stars in populations older 
than a few hundreds of Myr. For sake of consistency 
with the analysis in § 14.11 we consider only stars within 
1.5'. As already noted, even in this restricted area, field 
stars are distinguishable from cluster stars only in the 
youngest clusters, since in the remaining ones giant stars 
belonging to the cluster and to the field occupy the same 
regions of the CMD (see Figure [2]). Consequently, the 
contribution of field stars to the numerator can be eas- 
ily subtracted for the first six clusters. In older clusters 
we used the control-field frames, assuming that the field 
stellar population present in the adjacent field is identi- 
cal to that populating the cluster's area. The procedure 
was applied to each cluster of the sample and results are 
reported in Table [3] The errorbars also include the ef- 
fect due to the uncertainty on the membership of bright 
stars. 

We calculated theoretical SBFs from the exact simu- 
lations used to reproduce the CMDs, CLFs, integrated 
magnitudes and colors. In practice, for a given cluster 
SBF amplitudes were derived applying Eq. [5] and [5] In 
Figures [H and [TU] we plot differences between observed 
and predicted SBFs (model A and model C, respectively) 
as a function of the observed Mj and Mk s magnitudes. 
Differently from colors, that are reproduced by both as- 
sumptions (even if with different precision), the SBF 
data/models comparison shows a clear different behav- 
ior. SBFs from model A provide a good agreement with 
data in the entire range of Mj and Mk s , being differ- 
ences mostly lower or of the order of ~ 0.5 mag. There 
are only three outlayers, namely NGC 2134 and, at less 
extent, NGC 1783 and NGC 1806, in which a high field 
contamination was recognized (Figure [3]) . On the other 
hand, model C foresees SBFs fainter than the observed 
ones, in particular this model underestimates the number 
of bright stars in the magnitude intervals —8 ^ Mj ^ — 9 
and -9 < M K < - 9.8 (Figure Q35) . 

To better interpret these results, Figure [TT1 shows the 
relationship between SBF differences and s-parameter 
(age). Again, in the case of model A the agreement 
is good for all the s-parameter values considered here. 
On the other hand, it is clear that clusters with s > 30 
(ages greater than about 200 Myr), i.e. the age at which 
the number of TP- AGB stars becomes statistically rele- 
vant, suffer large discrepancy in the case of model C. Al- 
though the sample does not contain very old clusters, our 
guess is that for ages older tha n ~5 Gyr the effect shoul d 
decrease, as already shown in iRaimondo et all (|2005al) . 
Note that results from model B and D show a behav- 
ior similar to model C. Thus, changing the number of 
TP- AGB stars mainly affect the SWB classes IV, V and 



VI (see also Figure [5]). We recall the four models (.4- 
D) differ only for the treatment of the mass- loss process 
of TP- AGB stars. On this basis, we suggest that model 
A is more appropriate to reproduce the SBF measure- 
ments of LMC star clusters in the range 30 < s < 45, 
where TP- AGB stars are expected to have a large im- 
pact on the cluster observational properties. It is worth 
noting that this conclusion is not unique and should be 
verified when different assumptions are used in stellar 
population models (e.g. IMF, stellar evolutionary tracks, 
color-temperature relations, TP- AGB synthetic prescrip- 
tions, etc.). Within this limitation, we suggest that SBFs 
can contribute to disentangle the observable properties of 
TP- AGB stars in unresolved galaxies containing a large 
number of such stars. 

Before concluding this sectio n, we compare presen t 
SBF measurements to those by Gonzal ez et al.l (|2004f l. 
They derived NIR-SBFs of star clusters in the MCs as 
grouped (to o btain a sort of sup er-clusters) according to 
SWB classes (i Searle et al.| [l980). This procedure allows 
the authors to reduce the problem of small number statis- 
tics, and to compare LMC cluster data to SBF predic- 
tions computed for systems with a large number of stars, 
i.e. galaxies (see IRaimondo et aLll2005a! for details). We 
grouped the present clusters using the same definition of 
SWB classes. Our sample does not cover the Pre, SWB- 
I and VII classes. Moreover, the SWB-II class contains 
only one cluster (NGC2164), thus it is not significative. 
Table [4] summarizes the cluster grouping and the result- 
ing SBF data. Figure [12] compares SBF data of super- 
clusters from both sets of measures. Despite the smaller 
number of clusters considered in the present work, the 
agreement between the two sets of measurements is quite 
good, even though there is a systematic age shift of SWB- 
III and SWB-V due to t he diff erent s-parameter/age re- 
lation. [Gonzalez et al. (2004) assig ned ages and meta l- 
licities to the super-clusters from iFrogel et al.l (|1990l) . 
with a few exceptions regarding SWB-I, SWB-II and Pre- 
SWB super-cluster. The open stars in Figure [12] result 
from grouping the corresponding synthetic CMDs, i.e. 
we applied exactly the same procedure used for observa- 
tions. The proper SWB class is assigned according to the 
present age estimations of each cluster, then a mean is 
considered. We note a slight discrepancy (<0.5 mag) for 
SWB- VI in the J-band and for SWB-V in the J^-band. 
We recall that for these two classes the overlapping be- 
tween cluster and field population is more severe. As 
a conclusion, we note that the data/models comparison 
for super-clusters stands again for a consistency of the 
present procedure and SPS models. 

5. SUMMARY AND CONCLUSIONS 

In this paper we carried out a joint analysis of resolved 
and unresolved properties of a sample of LMC clusters, 
in order to consistently derive the evolutionary status of 
stars in the clusters, by using different stellar population 
synthesis tools based on an unique theoretical framework. 
The final aim was understanding the potential of the SBF 
technique as tracer of stellar populations, and exploring 
the use of such a technique to constrain more tightly 
stellar population models, those of intermediate age in 
particular. 

To reach this goal, we studied the CMD, CLF, in- 
tegrated magnitudes and colors, and SBFs of 19 LMC 
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star clusters, covering an age range from a few Myr 
up to several Gyr. We used accurate NIR (J and 
K s bands) photome t ric data of cluster's res olved stars 
(jFerraro et al.l 12004 iMucciarelli et all 120061), and NIR 
colors from up-to-date com pila tions by [G oudfrooii et alJ 
(p006T ). lPessev et al.l (f2006f ) and IMucciarelli et all (12006). 
By coupling resolved star photometry with total fluxes, 
we derived new measurements of NIR-SBFs for each indi- 
vidual cluster, paying attention to field star contamina- 
tion which may highly affect stellar counts in the upper 
part of the luminosity function. This allows us to obtain 
a set of homogeneous and high-quality SBF data of pure 
intermediate-age stellar populations. This is a necessary 
step toward having a larger sample of star clusters with 
accurate measurements of SBFs. 

We used the multi-purpose synthesis code SPoT, based 
on a unique and consistent theoretical framework, to pre- 
dict all the quantities quoted above. Great attention 
was payed to take into account stochastic effects. This 
is a crucial point, as broad-band colors and, in general, 
spectro-photometric features of clusters may suffer from 
large intrinsic fluctuations cau sed by the discrete nature 
of the number of s tars (see e.g.lBarbaro fc Bertellilll977l: 
Chiosi et all 119881: ISantos fc Frogel l 119971: iBro cato et all 
1999at iRaimondo et al.l l2005al) . Therefore, including 
stochastic effects is an important step to derive reliable 
age and m etallicity estimations from cluster broad-band 
colors (e.g. iFagiolini et aLll2007h . 

We also presented new synthesis models (A-D in Sect. 
12. 2|) . computed under the same prescriptions for all the 
evolutionary phases except for the mass-loss rates of 
TP-AGB stars. Model results were compared to data 
through the fitting of observed CMD, CLF, and inte- 
grated colors of each cluster. The agreement obtained 
in reproducing both resolved and unresolved observa- 
tional quantities suggests a high degree of reliability of 
the adopted theoretical framework. The cluster ages and 
metallicities derived by best-fit models are also in agree- 
ment (within the uncertainties) with previous determi- 
nations from other indicators. The age estimates were 
used to derive a new calibration of the s-parameter/age 
relationship for intermediate-age stellar clusters, i.e. s 
values within the range 20-45. 

Follo wing the same line presented in IRaimondo et al.l 
(2005a), we showed that NIR-SBFs are extremely sen- 
sitive to the number and properties of AGB (especially 
TP-AGB) stars. Changing the mass-loss rate prescrip- 
tion affects TP-AGB properties, including the duration 
of the entire phase and thus SBF predictions change siz- 
ably. On the basis of the theoretical framework adopted 
(WG98 and W90) to describe the TP-AGB evolution, 
we found that although all the models are able to re- 
produce observed CMDs, CLFs and integrated magni- 
tudes and colors, even though with different precision, 
not all the models reproduce the SBF measurements of 
LMC clusters. In particular, only model A provides a 
good fit of NIR-SBF in the age range corresponding to 
30 < s < 45, where TP-AGB stars are expected to be 
more numerous. Thus, we argue that the global theo- 
retical scenario adopted in model A appears to be able 
to reproduce all the observational properties of the LMC 
star clusters investigated in this work. It is worth notic- 
ing that the theoretical prescriptions used for modeling 



stellar clusters may affect our results on TP-AGB syn- 
thetic evolution, for instance the luminosity and effective 
temperature at the first thermal pulse. In this sense, an 
effort to study the physics of AGB structure experienc- 
ing thermal pulses is urged to provide solid theoretical 
background for stellar population synthesis studies. In 
spite of this limit, we showed that NIR-SBFs are very 
sensitive to the properties of TP-AGB stars, then they 
can be useful to investigate the global nature of such 
stars in clusters and galaxies, where they are numerous 
and statistical effects are less important. Using the large 
and homogeneous sample of clusters presented here, we 
have shown that SBF amplitudes may be linked to stellar 
evolution properties, and, in turn, they can be used as 
an intriguing and complementary diagnostic to test the 
ingredients of SPS models, and evolution properties of 
bright giant stars. 

As conclusive discussion, we notice that the metal- 
licity dependence of the C-rich phase lifetime was ac- 
counted for using the relationship between the number 
ratio of C- to M-stars and the metallicity obse rved for 
the Local Group galaxies (jGroenewege a l2007l ). How- 
ever, the number ratio of C-rich/O-rich stars should be 
a complex function of the properties of individual AGB 
stars, i.e. the durations of the entire phase and the sub- 
phases (O/C-rich, respectively), and the location of AGB 
stars in the CMD. Moreover, even though the use of a 
purely synthetic code to describe the TP-AGB phase 
is less accurate than the full stellar modeling, it is the 
only way to explore a large parameter space in com- 
puting SPS models. The uncertainties in assigning - 
as well as in predicting- stellar parameters of variable 
TP-AGB stars (e.g. effective temperature, opacity, lu- 
minosity), where effects such as convection, variability, 
mass los s, and dust formation b ecome increasingly im- 
portant (|Frevt ag & Hofncr 2008), still have a great im- 
pact in SPS models. Modern stellar evolutionary codes 
as well as hydrodynamical codes for stellar atmospheres 
are at present extensively used to model different aspects 
of stellar physics at work in stars beyond the core-helium 
burning phase, as the interior nucleosynthesis and enve- 
lope chemical enrichment and their dependence on stel- 
lar metallicity and mass, hence we expect to reduce the 
uncertainty of the predicted properties of AGB stars and 
then of intermediate-age stellar population models in the 
next future. Once again detailed analysis and accurate 
simulations of well-known stellar populations like LMC 
star clusters represent a fundamental step to improve 
our understanding of the nature of contributors to the 
integrated light from remote galaxies, as well as the po- 
tential of SBFs to disentangle the evolutionary status of 
unresolved populations. 
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Fig. 6. — Integ rated magnitudes and colors for 19 LMC clusters. Predicted valu es refer to Model A. Upper left panel: Predicted My vs. 
observed My of van den Bcrgh (1981, open squares), and Goudfrooii ct al. (2006, filled circles). Upper right panel: The Mk s magnitude 
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As in Figure [6] but for Model C. 
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Fig. 8. — A" 3 -magnitude and J — K s color differences for model Model B (upper panels) and Model D (lower panel). 
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Fig. 9. — The difference between observed and predicted SBF magnitudes. Symbol size decreases from SWB class II to VI as labeled. 
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Fig. 10.— 



As in Figure [9] but for Model G. 
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TABLE 1 
LMC Cluster Sample Information 



Cluster Name 


Log Age 


Reference 


rev, /xji 
|*e/±i| 


Reference 


C-stars 


Reference 


s 


(!) 


(ucx ) 
( 2 ) 


(3) 


(acx ) 
( 4 ) 


(5) 


(Number) 
(6) 


(7) 


(8) 


IN uU ±DO± 


q o n +0.08 
9.08 



1 


fl "37-1-0, 9(1 

— u.o / nzu.zu 

— U.oo±U.Uo 


7 

12 


1-1-1 

1,1,1 


1 -Q-1 1 

i ,y , 1 1 


ay 


JNCjL/ l/oo 


8.94 


1 


— 0.45 


8 


2;4;4 


1)9)11 


3 i 


NGC 1806 


9.16 


1 


-0.23 


7 


4;2 


1;9 


10 


JNtjO 1831 


o c n I n on 

8.50±0.30 
8.58 


1 
1 


r\ m in on 

0.01±0.20 


7 


3; 2 


1)11 


31 


IN^rO loOO 


q i o_i_n Qfl 

8.20 


1 

1 

1 


n e;n+n 1 n 


a 




n ■ n 
u.u 


1 .Q 


7 


NGC 1978 


9.52 


1 


-0.96 


6 


4;4;9 


1;9;11 


15 


NGC 1987 


8.79 


1 


-1.00 


8 


3;1;1 


1;9;11 


35 


NGC 2031 


8.20±0.10 
8.20 


2 
1 


-0.52±0.21 


2 





1 


27 


NGC 2108 


8.86 


1 


-1.20 


8 


1;1;1 


1;9;11 


36 


NGC 2134 


8.28 


1 


-1.00 


8 


0;0 


1;11 


28 


NGC 2136 


8.00±0.10 
8.13 


2 
1 


-0.55±0.23 


2 


1;0;0 


1;9;11 


26 


NGC 2157 


i .60±0.20 
8.06 


3 
1 


—0.45 
-0.60 


10 

8 





1 


25 


T\T /"< f~l 01 CO 

InCjC 2162 


o 11+0.12 
9-H-0.16 

9.08 


5 
1 


n no I n nn 

— 0.23±0.20 
-0.46±0.07 


7 
12 


1)0 


1)11 


39 


NGC 2164 


7.70±0.20 
7.91 


i 
1 


—0.45 
-0.60 


10 

8 





1 


23 


NGC 2173 


Q no+0.07 

9.30 


6 
1 


-0.24±0.20 
-0.42±0.03 


7 
12 


1;1;1 


1;9;11 


12 


NGC 2190 


8.86 


1 


-0.12 


7 


2:2 


l;ll 




NGC 2209 


8.98tg;5S 
8.79 


5 
1 


-0.47 
-1.20 


10 

8 


2;2;2 


1;9;11 


35 


NGC 2231 


9 18+° 10 
8.94 


5 
1 


-0.67±0.20 
-0.52±0.03 


7 
12 


1;1;1 


1;9;11 


37 


NGC 2249 


8.82±0.30 
8.72 


i 
1 


-0.47 
-0.12 


10 

8 





1 


34 



Note. — Refere nces: (1) IMucciaxelli et all (120061), (2) IDirsch et all <2000f). (3) lElsonl dl99]f) . (4) 
Elson fc Fafll U988I), ( 5) IGeisler et all (1997), (6) Hill ct a . (2 0~0C|), (7 ) IQlszewski et aU 1119911). (8) 
Sagar fc Pandevl HlQgSTI. C91 IFrogel et al.1 1 119901) . (10) IMackev &; Gilmorei (120031) . (11) IWesterlund et al.l 
H991B . (12) IGrocnolski et al.1 (120061) . 



TABLE 2 

Synthetic properties of cluster sample 



Cluster 


log Age 




Z 




My 




M Ks 




-K S 




H 


-K s 




Name 


(yr) 




(dcx) 


(mag) 




(mag) 


( 


mag) 




( 


mag) 




(1) 


(2) 




(3) 




(4) 




(5) 




(6) 






(7) 




NGC 2164 


7.70+.°, 


18 
15 


2.0e- 


02 


-8.42±0 


10 


-9.7±0.3 





49±0 


09 





12±0 


05 


NGC 2157 


7.72+.°, 


17 
15 


2.0e- 


02 


-8.68±0 


06 


-10.0±0.3 





52±0 


11 





13±0 


07 


NGC 2136 


7.72+.°, 


17 

37 


2.0e- 


02 


-8.62±0 


09 


-9.9±0.3 





50±0 


11 





14±0 


07 


NGC 1866 


8.08+.° 


07 
10 


8.0e- 


03 


-9.07±0 


03 


-10.8±0.3 





61±0 


15 





17±0 


09 


NGC 2031 


8.08±0 


OS 


2.0e- 


02 


-8.07±0 


05 


-9.9±0.6 





67±0 


28 





22±0 


16 


NGC 2134 


8.18l° 


05 
10 


2.0e- 


02 


-7.99±0 


04 


-10.2±0.6 





75±0 


25 





25±0 


15 


NGC 1831 


8.50±0 


07 


2.0e- 


02 


-8.07±0 


02 


-10.1±0.4 





79±0 


18 





28±0 


11 


NGC 2249 


8.68±0 


09 


8.0e- 


03 


-6.69±0 


06 


-8.7±0.6 





68±0 


26 





22±0 


11 


NGC 2190 


8.81±0 


07 


8.0e- 


03 


-6.73±0 


06 


-9.5±0.6 





93±0 


26 





33±0 


15 


NGC 2209 


8.95±0 


10 


8.0e- 


03 


-6.33±0 


07 


-9.2±0.7 


1 


01±0 


25 





34±0 


15 


NGC 1987 


8.95±0 


10 


8.0e- 


03 


-7.04±0 


05 


-10.1±0.4 


1 


08±0 


14 





37±0 


10 


NGC 2108 


8.95±0 


10 


8.0e- 


03 


-7.28±0 


05 


-10.1±0.5 


1 


06±0 


17 





37±0 


11 


NGC 2231 


8.98±0 


11 


8.0e- 


03 


-6.29±0 


08 


-9.3±0.6 


1 


02±0 


21 





34±0 


15 


NGC 1783 


9.08+.°, 


15 
23 


8.0e- 


03 


-8.28±0 


01 


-11.1±0.3 


1 


05±0 


13 





36±0 


10 


NGC 2162 


9.18±0 


Hi. 


8.0e- 


03 


-6.44±0 


06 


-9.2±0.5 





88±0 


22 





27±0 


11 


NGC 1651 


9.23+° 


07 
12 


8.0e- 


03 


-6.60±0 


06 


-9.6±0.4 





97±0 


15 





32±0 


11 


NGC 1806 


9.34+.° 


10 

14 


4.0e- 


03 


-7.72±0 


03 


-10.4±0.2 





83±0 


11 





23±0 


09 


NGC 2173 


9.48±0 


08 


8.0e- 


03 


-6.89±0 


03 


-9.7±0.3 





89±0 


15 





26±0 


11 


NGC 1978 


9.70±0 


08 


4.0e- 


03 


-8.95±0 


01 


-11.3±0.2 





65±0 


05 





10±0 


07 
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TABLE 3 

SBF MAGNITUDES OF LMC CLUSTERS 



Models Observations 



Cluster Name 


M T 


lvl K s 


M T 


lvl K s 






(mag; 


[ XY\ Q <T 1 

(jnag; 


( VY1 Q(T 1 

(jnag; 


^mag; 


(1) 




V>! 


(d\ 
W 


C5~l 


NGC 


2164 


-5.1 ±0.6 


-6.1 ±0.8 


-5.1 ± 0.2 


-5.7 ± 0.8 


NGC 


2157 


-5.3 ± 1.0 


-6.4 ±1.3 


-5.6 ± 1.4 


-6.5 ± 1.3 


NGC 


2136 


-5.2 ±0.7 


-6.3 ± 1.0 


-5.0 ± 1.5 


-5.1 ± 1.6 


NGC 


1866 


-5.5 ±0.9 


-7.0 ± 1.2 


-5.5±0.1 


-7.0 ±0.1 


NGC 


2031 


-5.5 ± 1.8 


-7.0 ±2.2 


-6.4 ±0.1 


-7.6 ±0.3 


NGC 


2134 


-5.8 ± 1.6 


-7.4 ± 1.9 


-8.5 ±0.1 


-10.3 ± 0.1 


NGC 


1831 


-5.5 ±1.0 


-7.4 ± 1.2 


-4.8 ±0.3 


-6.9 ±0.1 


NGC 


2249 


-4.5 ± 1.8 


-5.9 ±2.2 


-4.6 ± 1.5 


-6.0 ± 1.0 


NGC 


2190 


-5.7 ± 1.0 


-7.5 ± 1.3 


-5.7±0.1 


-7.8 ±0.1 


NGC 


2209 


-5.6 ± 1.4 


-7.4 ± 1.7 


-5.6 ±0.2 


-7.6 ±0.1 


NGC 


1987 


-5.9 ±0.5 


-7.9 ±0.6 


-6.1 ±0.9 


-7.8 ±0.6 


NGC 


2108 


-5.9 ±0.6 


-7.8 ±0.7 


-6.2 ±0.4 


-8.0 ±0.1 


NGC 


2231 


-5.6 ± 1.0 


-7.4 ± 1.3 


-6.7 ±0.8 


-7.5 ±0.2 


NGC 


1783 


-5.6 ±0.5 


-7.7 ±0.5 


-4.7 ± 0.1 


-6.3 ±0.1 


NGC 


2162 


-5.0 ± 1.0 


-6.7 ±1.4 


-6.0 ±0.4 


-8.3 ±0.6 


NGC 


1651 


-5.3 ±0.7 


-7.2 ±0.9 


-5.7 ± 0.3 


-6.9 ±1.1 


NGC 


1806 


-5.0 ±0.5 


-6.9 ±0.6 


-6.9 ±0.1 


-7.1 ± 0.1 


NGC 


2173 


-4.7 ±0.7 


-6.6 ± 1.0 


-5.2 ±0.1 


-6.5 ±0.1 


NGC 


1978 


-4.3 ±0.1 


-5.4 ±0.1 


-4.7 ±0.1 


-6.4 ±0.4 



TABLE 4 

Cluster grouping and SBF magnitudes 









Models 


Observations 


SWB Class 


s 


Cluster Name 


Mj 




Mj 




II 


22-24 


2164 










III 


25-29 


2157,2136,1866,2031,2134 


-5.8 


-7.5 


-5.7 ± 0.7 


-7.1 ±0.4 


IV 


30-36 


1831,2249,1987,2209,2190,2108 


-5.6 


-7.7 


-5.6 ±0.5 


-7.6 ±0.2 


V 


37-41 


2231,1806,1651,2162,1783 


-5.4 


-7.5 


-5.4 ±0.6 


-7.0 ±0.3 


VI 


42-46 


2173,1978 


-4.4 


-6.1 


-4.7 ±0.2 


-6.4 ±0.3 



